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Abstract 



The specific heat of the x — y modei is studied on cubic iattices of sizes L x L x L and on lattices 
L x L x H with L ^ H (i.e. on lattices representing a film geometry) using the Cluster Monte Carlo 
method. Periodic boundary conditions were applied in all directions. In the cubic case we obtained 
the ratio of the critical exponents ajv from the size dependence of the energy density at the critical 
temperature T\. Using finite-size scaling theory, we find that while for both geometries our results scale 
to universal functions, these functions differ for the different geometries. We compare our findings to 
experimental results and results of renormalization group calculations. 

1 Introduction 

Physical systems which exhibit a second order phase transition and are confined in a finite geometry (e.g. 
a cubic or film geometry) are thought to be well described by the finite-size scaling theory at temperatures 
close to the critical temperature T\ [1]. This scaling theory states that finite-size effects can be observed 
when the bulk correlation length £ becomes of the order of the finite extent L of the system, e.g. in a cubic 
geometry the length of the edges of the cube plays the role of L. For a physical observable O this intuitive 
assumption can be cast into the following formula [2]: 



where £(t, L = oo) is the correlation length for the infinite size system, t is the reduced temperature, and / 
is a universal function. For example, the most singular behavior of the correlation length close to the critical 



where the prefactor £ has been absorbed in the definition of the universal function f^(x). 

A physical system, which has been widely used to experimentally test the finite-size scaling theory, is 
liquid 4 He since the superfluid density p s and the specific heat c can be measured to a very high accuracy. 
However, the experimental verification of the finite-size scaling theory is somewhat controversial. Rhee, 
Gasparini, and Bishop measured the superfluid density of very thick helium films [3] and showed that the 
data did not follow the form (1) using a value of v reasonably close to the expected value v = 0.67. Similarly, 
in early measurements of the specific heat of 4 He in finite geometries very different critical exponents from 
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the expected values were found [4] indicating that phenomenological finite-size scaling cannot be applied in 
a straightforward way. In order to clarify the situation renormalization group calculations for the standard 
Landau-Ginzburg functional in different geometries with Dirichlet boundary conditions have been under- 
taken [5, 6, 7, 8]. New specific heat measurements [9] and also a reanalysis [10] of the old specific heat data 
[4] show good agreement between the renormalization group calculations reported in [5, 6, 7] and those data. 
Furthermore new experiments on liquid 4 He confined in a film geometry under microgravity conditions are 
planned [11] in order to determine the finite-size scaling of the specific heat. Thus, it is desirable to obtain 
reliable results for the specific heat in finite-size helium systems by means of numerical investigations. 

In this paper we perform a numerical study of the scaling behavior of the specific heat of 4 if e in a cubic 
(L x L x L size lattices) and in a film geometry (L x L x H size lattices with L >> H) at temperatures close 
to the critical temperature T\. Since 4 He belongs to the universality class of the x — y model [12], we use 
this model to compute the specific heat at temperatures near T\ using the 1-cluster Monte Carlo method 
[13]. The x — y model on a lattice is defined as 



where the summation is over all nearest neighbors, s = (cos sin #) is a two-component vector constrained 
to the unit circle, and J sets the energy scale. 

The critical exponents of the three-dimensional x — y model have been determined by high-temperature 
expansions [14] and Monte Carlo simulations [15, 16, 17, 18]. The importance of vortex lines for the phase 
transition was investigated in Ref. [19]. A renormalization group approach based on vortex lines [20] 
derives the critical properties of the three-dimensional x — y model from the interaction of vortex lines. 
The anisotropic three-dimensional x — y model (J x = J y ^ J z ) has also been studied [21], and a crossover 
from three-dimensional to two-dimensional behavior was found with respect to the ratio J z / J x . The Villain 
model, which is in the same universality class as the x — y model, has been studied in a film geometry 
where the correlation length in the disordered phase was used to extract the thickness-dependent critical 
temperature [22]. The authors of Ref . [23] computed the universal scaling function of the superfluid density 
of helium confined in a film geometry using the x — y model and examined the crossover properties from 
three to two dimensions of the superfluid density. 

In this work we compute the specific heat c(T, L) of the x — y model on various cubic lattices L x L x L. 
We deduce its critical exponent a from the size-dependence of the energy density at the critical temperature, 
estimate the bulk value c(7\,oo), and check the scaling hypothesis for the specific heat with respect to L. 
We compare the resulting universal function to recent renormalization group calculations [26]. The specific 
heat for a film geometry, i.e. for various L 2 x H lattices with L — ► oo, is also computed and the scaling 
behavior of the specific heat with respect to H is studied. We compare the universal scaling function for the 
film geometry directly to the experimental results of Refs. [9]. 

The article is organized as follows. In the next section we introduce the definition of the specific heat and 
the energy density and briefly discuss the Monte Carlo method. Section 3 discusses the finite-size scaling 
properties of the specific heat. In section 4 we deduce the critical exponents and check the scaling assumption 
of the specific heat for the cubic geometry. Section 5 is devoted to the film geometry, and the last section 
summarizes our results. 
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2 Definition of the physical quantities and Monte Carlo method 



We define the energy density of our model as follows: 



E = (e) = 3 - 




where V = L 3 for the cubes and V = H L 2 for the film geometry. The specific heat can be written as 
(k B = 1): 

c=Vf({e 2 )-{e) 2 ), (5) 

where [3 = J /{hsT). 

The thermal averages, denoted by the angular brackets, are computed according to 

i 

0[6] denotes the dependence of the physical observable O on the configuration {0i}, and the partition function 
Z is given by 

^ = /n^exp(--^). (7) 

The multi-dimensional integrals in the expressions (6) and (7) are computed with Wolff's 1-cluster algo- 
rithm Monte Carlo method [13]. We computed the specific heat at various temperatures on _L 3 -lattices for 
L = 20, 30, 40 and on L 2 x if -lattices for L = 40, 60, 100 and H = 6,8, 10. At the critical temperature the en- 
ergy density and the specific heat were computed on L 3 lattices for L = 10, 15, 20, 25, 30, 35, 40, 45, 50, 60, 80. 
Periodic boundary conditions were applied in all directions. We carried out of the order of 20,000 thermal- 
ization steps and of the order of 500,000 measurements. The calculations were performed on a heterogeneous 
environment of computers including Sun, IBM RS/6000 and DEC alpha AXP workstations and a Cray-YMP. 



3 Finite— size scaling properties of the specific heat 

Let us first consider the x — y model in a cube whose edges are of length L. In such a geometry the specific 
heat starts feeling the finite size of the cube when the bulk correlation length £ becomes comparable to the 
length L. At temperatures close to the critical temperature T\ and for L — ► oo we can write [2] 

44 = (8) 

c(t, oo) 

where the reduced temperature t = T/T\ — 1, G is a universal function, and v is the critical exponent of the 
correlation length. In the limit L — ► oo and at a fixed value oft we obtain G(±oo) = 1. Now we leave L 
fixed but carry out the limit t — ► ±0 assuming c(t, oo) oc \t\~ a with a > 0. We obtain 

t lim c(t,L)oc t lim Itl-^L 1 /"), (9) 

and the fact that c(0, L) is finite at t = implies that 

lim G{x) cx \x\ a , (10) 

x— ^±0 
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thus 

c(0, L) oc L a l" . (11) 

Experiments on superfluid 4 He [24, 25] indicate that v = 0.6705, and via the hyperscaling relation a = 
2 — 2>v = —0.0115 < 0, thus c(0, oo) is finite. In order to write c(t, L) in a scaling form similar to (8), notice: 

c(t, oo) - c(0,oo) oc \t\~ a , (12) 

with a < 0. Since the scaling theory deals only with the most "singular" terms of a physical quantity when 
the critical point is approached, the following scaling form suggests itself: 

C (f,£)- C (0,oo) = lf 

c(t, oo) - c(0,oo) K J K ' 

or 

c(t,L) = c(0,co)+\t\- a g(tL 1 ^). (14) 

Keeping t fixed, we find lim^^-too g(x) = cf , which is finite. For fixed L and t — ► ±0 we obtain the behavior 
of Eq. (10) and 

c(0,L) = c(0,oo) + Cl L a/v . (15) 
After defining the scaling function g(x) = |a;| _a (jr(a;) we obtain 

c(t,L) = c(0,oo) + L a ^g(tL 1 ^). (16) 
This enables us to reexpress the scaling form (13) as follows: 

c(t,L) - c(0,oo) i; 
c(0,L)-c(0,oo) " Gi(<L } - 

Note that G L (0) = 1. 

Let us finally derive the relationship between the scaling function G^th 1 ! 11 ) and the scaling function 
fiitL 1 /") used in Refs. [5, 6]. This function is defined as follows: 

c(t,L)-c(to,^) = L a ^f 1 (tL 1 ^), (18) 

where t = (^ /L) 1 !" , i.e. t is the reduced temperature where the correlation length is equal to the system 
size L. Using Eqs. (15) and (17) we can write 

c(t, L) - c(0, oo) = c 1 L a l v G L {tL 1 l v ). (19) 

In order to make Eq. (19) consistent with Eq. (10), we have to require that if we keep t fixed and take the 
limit L — ► oo 

lrni G L (x) = gt\x\- a , (20) 

where g^ is a constant. Thus, at L = oo Eq. (19) can be written as: 

c(t,oo)-c(0,oo) = cf\t\- a , (21) 

where 2^ = cig^. Evaluating expression (21) at to, solving for c(0, oo) and inserting the result into Eq. (19) 
yields 

c(t, L) - c(t , oo) = L a ' v Cl {GdtL 1 '") - gUti)-" 1 "} ■ (22) 
Comparing this expression and the definition for fi(x) given by Eq. (18) we find 

fox) = Cl {G L (x)-g+ (£+)-"/"}. (23) 

These scaling forms are valid for the film confining geometry also. In the case of a film geometry we need 
to replace L by H in the scaling equations (13), (17), and (23) because the relevant scale is H when L H . 



4 



4 The cubic geometry 

In this section we investigate the finite-size scaling behavior of the specific heat of the x — y model in the 
cubic geometry L 3 . 



-0.5 -0.3 -0.1 0.1 0.3 0.5 

10t 



Figure 1: The specific heat for various size lattices I as a function of the reduced temperature. The solid 
lines represent L = oo deduced from our Monte Carlo calculation. T\ = 2.2017. 

In Fig.l we show our data for the specific heat. As a comparison we also plotted the bulk behavior of 
the specific heat (solid line). The steps leading to this curve are given below in this section. 



L 


c(0,L) 


E(0,L) 


10 


1.9454(69) 


1.95435(32) 


15 


2.180(11) 


1.98078(22) 


20 


2.362(10) 


1.99176(18) 


25 


2.465(13) 


1.99709(21) 


30 


2.579(18) 


2.00059(13) 


35 


2.621(19) 


2.00257(15) 


40 


2.754(30) 


2.00464(28) 


45 


2.789(34) 


2.00522(13) 


50 


2.783(31) 


2.00618(10) 


60 


2.967(45) 


2.00712(10) 


80 


3.048(43) 


2.00882(11) 



Table 1: The Monte Carlo results for the specific heat c(0, L) and the energy density £"(0, L) at the critical 
temperature T\ = 2.2017. 

In order to find the universal function Gl{x) we need to know c(0,oo). This quantity can be found by 
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calculating the specific heat at the critical temperature for various lattice sizes L and fitting the data to the 
form (15). We take T\ = 2.2017 as the critical temperature [18]. Table 1 contains the values for c(0, L) and 
for £"(0, L). Since the specific heat is very sensitive to fluctuations, it has a relatively large error, making it 
very difficult to extract the very small exponent ajv . Therefore we decided to use another procedure to find 
the values of the critical exponents. A quantity which is closely related to the specific heat is the energy 
density E. It is advantageous to use the energy density data because of the small error bars involved in its 
calculation. These error bars are two orders of magnitude smaller than the error bars of the specific heat. 
From the expression 

«.L ) = »J%i>, (24 , 

we obtain by integrating (16) up to a constant 

E(t, L) = c(0, oo)T + L (a - 1 ^ v T x D{tL 1 l v ), (25) 
where dD(x)/dx = g(x). For t — ► we obtain 

E(0,E) = E + E 1 L (a ~ 1 V v . (26) 



data points 


E 


E 1 




ajv 


x 2 


Q 


10 


2.0111(2) 


-1.80(15) 


1.479(40) 


-0.029(20) 


1.58 


0.15 


9 


2.0110(4) 


-1.81(38) 


1.487(81) 


-0.0258(75) 


1.76 


0.12 


8 


2.0111(7) 


-1.80(10) 


1.48(25) 


-0.029(55) 


1.92 


0.10 


7 


2.011(2) 


-1.8(37) 


1.48(76) 


-0.029(88) 


2.56 


0.05 



Table 2: Fitted values of the parameters entering expression (26). Q is the goodness of the fit. 



data points 


c(0, oo) 


Cl 


x 2 


Q 


10 


21.33(50) 


-20.57(54) 


1.51 


0.15 


9 


20.45(66) 


-19.61(72) 


1.12 


0.35 


8 


20.72(94) 


-19.9(10) 


1.28 


0.26 


7 


20.2(13) 


-19.3(15) 


1.47 


0.20 



Table 3: Fitted values of the parameters entering expression (15). ajv = —0.0258. Q is the goodness of the 



The results of the fits of the energy density data to the expression (26) are given in Table 2. In the fits 
we have subsequently excluded values of the energy density corresponding to smaller and smaller lattices. 
Because of the size of the error bars of the ratio ajv and 1/v we cannot make a definite statement that 
ajv < and \jv < 1.5. The fitting parameters, however, become stable when the energy density data 
obtained for lattice sizes L > 20 (9 data points) are used for the fits. Despite the large error bars, we always 
find ajv < and \jv < 1.5. The parameters for the fit including 9 data points and shown in Fig. 2 are 

E = 2.0110 ± 0.0004, 

E l = -1.81 ±0.38, 

\jv = 1.487 ± 0.081, (27) 

ajv = -0.0258 ± 0.0075. (28) 
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Figure 2: The energy density £"(0, L) at the critical temperature T\ = 2.2017 as a function of L. The solid 
curve represents the fit to (26) (9 data points included). 
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Figure 3: The specific heat c(0,L) at the critical temperature T\ = 2.2017 as a function of L. The solid 
curve represents the fit to (15) (9 data points included), ajv = —0.0258. 
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Within error bars the hyperscaling assumption is fulfilled. The experimental value for ajv is ajv = —0.0172 
[24], and an earlier experiment gave ajv = —0.0225 [25]. Having determined ct/v, we can turn to fitting 
the specific heat data to the expression (15). We fixed the value of ajv to the previously determined value 
ajv = —0.0258. Table 3 contains the fitting results. If we exclude the data corresponding to the two smallest 
lattices we obtain 

c(0,oo) = 20.45 ± 0.66, (29) 
ci = -19.61 ± 0.72, (30) 

and the fit is shown in Fig. 3. 




Figure 4: The scaling function — Gl{x) (cf. Eq. (17)) for the cubic geometry. c(0, oo) = 20.45, \jv = 1.487, 
T\ = 2.2017, the values for c(0, L) are taken from Table 1. 

In order to check the validity of the scaling form (17) we plot (c(t, L) — c(0, oo))/(c(0, oo) — c(0, L)) versus 
tL 1 !" for different lattice sizes L 3 in Fig. 4. We used the values for c(0,L) given in Table 1 and c(0,oo) as 
determined above (Eq. (29)). As expected the data points for the three lattices 20 3 , 30 3 , and 40 3 collapse 
onto one universal curve Gl^L 1 !") in the range -10 < tL 1 !" < 10. 

It is interesting to repeat the fits described above using the experimentally determined critical exponents 
v = 0.6705 and a = —0.0115 [24]. The result of the fit of the specific heat data corresponding to lattices of 
size L > 20 to the expression (15) is 

c(0,oo) = 30.3 ±1.0, (31) 
ci = -29.4 ±1.1. (32) 

(33) 

Fig. 5 shows the scaling plot where v = 0.6705 and c(0,oo) = 30.3. The data of the specific heat for the 
20 3 , 30 3 , and 40 3 lattice collapse onto one universal curve. We see, that the value of c(0,oo) is strongly 
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Figure 5: The scaling function — Gl{x) (cf. Eq. (17)) for the cubic geometry. c(0,oo) = 30.3, v = 0.6705, 
T\ = 2.2017, the values for c(0, L) are taken from Table I. 

effected by the value of the ratio ajv. A much wider range of lattice sizes is necessary in order to determine 
these values more accurately. However, the scaling property, i.e. the data collapse onto one universal curve, 
is rather insensitive to the precise value of the ratio ajv and thus to the precise value of c(0, oo). This is 
demonstrated in Figs. 4 and 5. Thus, for our lattice sizes the data collapse cannot be used to determine the 
value of the critical exponents more accurately. 

In Fig. 6 we compare the results of the renormalization group calculations of Ref. [26] (solid line) to the 
results of our Monte Carlo simulation for the lattice with L = 48. We computed the values for the specific 
heat on this lattice from our scaling function Gl{x) (cf. Fig. 5). The agreement between the renormalization 
group calculations and the Monte Carlo results is satisfactory. 

We now compute the function fi(x) defined by Eq. (23) for the x — y model. This function can be 
obtained from Eq. (23) with v = 0.6705, c\ = —29.4, and = 0.498 from Ref. [16]. We estimate the value 
for g+ as follows. Eq. (20) implies that if we plot x a GL{x) versus x for large enough positive values of x 
the function x a GL{x) should approach the finite value g+ . We obtain g+ = 1.0378(5) and thus 

h(x) = -29.4 [G L (x)- 1.025]. (34) 

This function is shown in Fig. 7. The solid line in Fig. 7 represents the result of the renormalization group 
calculation of Ref. [26] (the scaling function P c (x) given there is related to fi(x) by the relation fi(x) = 
Pc(%) + 19.7). Also here the agreement between the renormalization group calculation and the Monte Carlo 
results is satisfactory. We would like to note that the shape of the function fi(x) is rather insensitive to the 
precise value of ajv and thus to the precise value of c\ (as was the scaling property of the specific heat), i.e. 
the scaling function Gl{x) given in Fig. 4 yields almost the same function fi(x) as shown in Fig. 7. Therefore 
the scaling function fi(x) cannot be used to determine the critical exponents more accurately. However, the 
function fi(x) lends itself well to compare experimental results to the results of our calculation because in 
both cases it is hard to determine a and thus c\ to a better accuracy. 
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Figure 6: Comparison of the results of this work (MC) and the renormalization group (RG) calculations of 
Ref. [26] (solid line) for the specific on a L = 48 lattice. 
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Figure 7: The function fi(x) for the cubic geometry. The solid line is the result of the renormalization group 
(RG) calculation of Ref . [26], the squares represent our Monte Carlo (MC) simulation. 
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Figure 8: The specific heat for various lattice sizes L derived from the scaling function Gl{x)- The solid line 
represents L = oo according to Eq. (21). 
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Figure 9: The bulk specific heat as a function of the reduced temperature. The solid lines represent the 
results of this work, the dashed lines represent the results of high temperature expansions and of Ref . [26]. 
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In what follows we will determine the temperature dependence of the bulk specific heat above and below 
T\ using the experimentally determined values of v = 0.6705 [24] and a = —0.0115 and thus the values for 
the quantities c(0, oo) and c\ given by Eqs. (31) and (32). The knowledge of g+ enables us to find the bulk 
behavior of the specific heat above T\ as given by Eq. (21) where 

c+ = 015+ = -30.5 ± 1.1. (35) 

We found that below T\ the asymptotic form (21) is only valid for a narrow region near T\ (namely \t\ <C 
0.02). In Fig. 8 we show curves for the specific heat for lattice sizes ranging between L = 48 and L = 240. 
We computed these curves from our scaling function Gl{x) where a = —0.0115. The solid line in Fig. 8 is 
the bulk curve. We identified the bulk values of the specific heat by the values of the specific heat for large 
enough size lattices which collapse onto one curve, because at those temperatures the specific heat does 
not feel the finite size of such large systems (because the correlation length for those temperatures is much 
smaller than the size of these lattices). Using larger and larger lattice sizes we can reach temperatures closer 
and closer to T\. Very close to T\ the bulk curve can be expressed by Eq. (21) giving 

27/ = -29.225 ± 0.025, (36) 

where c(0,oo) is taken from Eq. (31). In Fig. 9 the bulk behavior of the specific heat is shown. The solid 
lines represent the results of the Monte Carlo calculation reported in this work. The dashed lines are taken 
from Ref. [26] which are a combination of the high-temperature series expansion of Ferer et al. (cf. Refs. 
[14]) and the renormalization group calculations of Ref. [26]. We would like to emphasize that the expression 
(21) is only valid in the interval -0.02 < t < 0. 
We can use the universal ratio [27, 28] 



c 



(37) 



to compare our results to experimental results. The experimental value of £|7 is £|7 = 3.57A where £|7 is 
defined from T/T(t) = £|7 \t\~" according to Refs. [28, 29] (T denotes the helicity modulus). Using the x — y 
model and the method of Ref. [23] we have calculated the helicity modulus for cubic lattices up to 40 3 and 
using finite-size scaling for T and the above definition of £|7 we find £|7 = 1.21 in lattice spacing units and 
thus a = 2.95A. With £+ = 0.498a from Ref. [16] we find £+ = 1.47A (£+ cannot be measured directly in 
the experiments). Since now the ratio £o"/£()" * s ^ ne same on the lattice and in the experiments so has to be 
the ratio 2^ /cj" . We have 

c+ _ f 1.058(4) from Ref. [30], 



Cl L 1.044(38) from this work. ^ 

The agreement is quite satisfactory though it would be desirable to determine this ratio more accurately on 
the lattice, which requires larger lattices and probably longer simulation times as were used in the calculations 
reported here. We can also compare = acf with the experimental results of Ref. [31] by expressing the 
specific heat of the x — y model in physical units. This is accomplished by the equation 

Joule 

Cs ~ a 3 C_14 -' °Kmole C ' 



where V m is the molar volume of 4 He at saturated vapor pressure at T\, c is given by expression (5), and 
c s has i 
obtain 



c s has units of Joule(°K mole) 1 . Using a = —0.0115 and the values for 2^ given by Eqs. (35) and (36) we 



A+ = 5.2 ± 0.02, (40) 
A~ = 4.954 ± 0.004, (41) 
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with in units of Joule(°K mole) 1 . The experimental results of Ref . [31] 



Also here the agreement is acceptable. 



A+ 
A~ 



-0.009, 

5.82, 

5.504. 



(42) 
(43) 
(44) 



5 The film geometry 

The specific heat on L 2 x H lattices in the limit L — ► oo should obey: 



c(t,H) - c(0,oo) 



c(0,#) - c(0,oo) 

The value of c(0, oo) is given either by (29) or (31) depending on the value for v. 



(45) 



2.5 



2.0 



= 1.5 



1.0 
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1.8 
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2.2 



2.4 



2.6 



Figure 10: The specific heat c(L, H , T) as a function of T for a film geometry for H 
of L. 



4 and different values 



In order to apply the scaling form (45) we need to take the limit L — ► oo. However we find that the 
Independence of the specific heat is very weak, already for L = 60 and L = 100 for a fixed H the specific 
heat data agree within error bars, as demonstrated in Fig. 10. This weak Independence of the specific heat is 
quite in contrast to the very strong L-dependence of the helicity modulus, whose values had to be computed 
in the limit L — ► oo at temperatures close to the U-dependent critical temperature T 2D [23]. We take 
the values of the specific heat computed on 100 2 x U-lattices for U = 6, 8, 10 and assume those values to 
represent the case of infinite planar dimension. 
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Figure 11: The scaling function —Gh(%) (cf. Eq. (45)) for the film geometry. c(0,oo) = 20.45, ajv 
—0.0258, 1/v = 1.487, T\ = 2.2017. The error bars have been omitted for clarity. 
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Figure 12: The scaling function —Gh(%) (cf. Eq. (45)) for the film geometry. c(0, oo) = 30.3, ajv 
v = 0.6705, T\ = 2.2017. The error bars have been omitted for clarity. 



-0.0172, 



14 



In order to verify the expression (45), we plot —Gh(%) in Fig. 11 where the values for the parameters 
a, v, c(0,oo), and c\ are taken from Eqs. (27), (28), (29) and (30). The data collapse onto one universal 
curve. Instead if we use the experimentally found critical exponents v = 0.6705 and a = —0.0115 and the 
values c(0, oo) and c\ given by (31) and (32), respectively, to check the validity of scaling we obtain also data 
collapse as shown in Fig. 12. The data points deviate from the scaling curve for x < —1.5 because we reach 
temperatures in this region which are outside the critical region. If we included data from much thicker films 
this deviation would occur for smaller values of x. 
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o D=H: film geometry 
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Figure 13: —Gh(%) for the film geometry (circles) and 
a/v = -0.0115, v = 0.6705, and c(0, oo) = 30.3. 



-Gl{x) for the cubic geometry (filled circles). 



Fig. 13 compares the universal functions Gl{x) for the cubic geometry and Gh(%) for the film geometry 
with v = 0.6705 and a = —0.0115. The universal function —Gn{x) has a sharper peak at its maximum. 

Finally in this section we would like to compute the function fi(x) in physical units for the film geometry 
in order to make a direct comparison with experiments. We determined the lattice spacing in physical 
units a and the value for g+ in the previous section is independent of the geometry and the boundary 
conditions). Thus, the function fi(x) in physical units (open circles) is determined up to a constant prefactor 
in front of the argument x which is fixed by requiring that the maximum of fi(x m ) occurs at the same value 
of the argument x m as the maximum of the experimentally determined function fi(x m ). This constant 
prefactor turns out to be 1 and fi(x) is displayed in Fig. 14. Also in this figure we plotted the experimentally 
determined function fi(x) of Refs. [9] (filled circles). The two functions are of the same order of magnitude 
but do not agree (varying the value of the constant prefactor does not improve the agreement). This is not 
surprising because the periodic boundary conditions only approximate the true physical boundary conditions. 
Therefore we also plotted the function fi(x) (solid line) obtained from a renormalization group calculation 
for the Landau-Ginzburg functional in a film geometry with Dirichlet boundary conditions [5, 6]. The 
importance of the correct physical boundary conditions is clear from this discussion. 
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Figure 14: Comparison of the experimentally determined function fi(x) (filled circles) of Refs. [9], fi(x) ob- 
tained for films with periodic boundary conditions (open circles), and fi(x) for films with Dirichlet boundary 
conditions of Ref. [6]. fi(x) has the units Joule(°K mole) -1 and H is in A. 



6 Summary 

We have investigated the finite-size scaling properties of the specific heat of the x — y model in a cubic 
geometry L x L x L and in a film geometry L x L x H with L H . Periodic boundary conditions were 
applied in all directions. For the cubic geometry we find strong evidence that the critical exponent a is 
negative, so the specific heat does not diverge at the critical temperature, in qualitative agreement with 
experimental findings. However, we were not able to determine a very accurately, we find ajv = —0.026(8), 
1/v = 1.49(8), which is in reasonable agreement with experiments. Our values for a and v fulfill the 
hyperscaling relation. We confirmed the scaling assumptions (17) for the cubic case and (45) for the film 
geometry. We also used the experimentally determined values for v and a to compute the scaling functions 
for the specific heat in the two geometries. We derived the bulk behavior of the specific heat below and 
above the critical temperature and compared these results to renormalization group calculations and high- 
temperature expansions. Good agreement between the scaling function fi(x) (cf. Eq. (18)) for the cubic 
geometry obtained from our Monte Carlo simulation and renormalization group calculations was found. For 
the film geometry we compared fi(x) derived from our Monte Carlo data to the experimentally determined 
function fi(x) and the renormalization group result of Refs. [5, 6]. This comparison leads to the conclusion 
that the boundary conditions determine the shape of the universal functions and have to be chosen properly 
in order to find agreement with the experimental results. 
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